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Abstract 

Finite-bandwidth effects on the temperature dependence of the mobility of injected 
carriers in pure organic crystals are explored for a simplifed case of impurity scat- 
tering. Temperature-dependent bandwidth effects are discussed briefly through a 
simplified combination of band and polaronic concepts. 
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1 Introduction 

In studying charge carrier (electron or hole) transport in organic crystals, 
one avoids the use of standard band theoretic approaches [1] since the carrier 
bandwidth is not an overwhelmingly large energy in organic, in contrast to 
inorganic, materials [2]. Mean free paths calculated on the basis of bare band- 
theoretic descriptions have been explicitly found to be smaller than a lattice 
constant [3]. Such findings have necessitated polaronic or hopping theories of 
transport in organic materials [4,5]. The question arises whether this trend 
[6,7] away from bare-band theory, which has continued in the literature on 
organic materials in recent times [8], needs to be revised or at least modified 
in light of recent observations [9] and calculations [10] on pentacene. The 
large low-temperature mobility in these materials (orders of magnitude larger 
than values reported earlier [11] at other temperatures and in other organics), 
sharp power-law dependence of the mobility on temperature in a wide range 
[9,12], and nonlinear saturation phenomena [13] which have been interpreted 
[9] in terms of large bandwidths, all could point to the need for bare-band 
theory to be considered seriously in organic materials. On the other hand, 
the bandwidths are believed [9,10] not to be of the order of several eV's as in 
inorganic materials but smaller. One concludes, therefore, that required is a 
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description valid for intermediate bandwidths (B). To begin the construction 
of such a description is the main purpose of the present Note. 

Temperature-dependent (Huang-Rhys) 5's arise in polaronic transport [4,5,14]. 
In the past, effective bare S's, which undergo reduction by polaronic effects in 
the presence of strong interactions with vibrations, were themselves believed 
(or shown) to be rather small. With naphthalene as an example, B ~ 10 eV, 
the polaronically reduced B is about a tenth of this value in light of the mag- 
nitude of the coupling constants and is thus well surpassed by the thermal 
energy k B T . The latter exceeds 2.5 meV corresponding to 30 K, the lowest 
temperature in the naphthalene observations [11]. Therefore, although past 
explanations of mobility observations have certainly made use of the tempera- 
ture dependence of polaronic bandwidths [4,5] to address decreasing mobility 
with increasing temperature [5], polaronic calculations have served as the in- 
vestigational procedure. Band considerations have either not been used at all 
[5] or assigned a less important role [4] . Estimates in pentacene raise the pos- 
sibility that k B T is smaller than B for a part of the observations. To show 
the results of a simple combination capable of blending finite-band theory and 
polaronic concepts, is a secondary purpose of the present Note. 



2 Constant Scattering Rates and Truncated Parabolic Bands 



Band-theoretic descriptions in inorganic materials effectively take the band- 
width B to be infinite. The well-known textbook formula for the mobility 

^ \k B Tj JS°dep(e)e-^BT W 

where k B is the Boltzmann constant, and q, e, T, v, p, r are the carrier 
charge, energy, temperature, velocity, density of states, and relaxation time, 
respectively, has oo as the upper limit of the energy integrations. We will 
explore, instead, the consequences of replacing it by the finite value of B, 
suitably modifying the density of states. We consider a one-band model, which 
means merely that the band gap is taken to be large enough to make interband 
transitions unimportant, and assume a simple p (e) having the free-electronic 
form within the band and to vanish outside: 

p{e)=N(^0=\[e{e)-9{e-B)} (2) 
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where 9 is the Heaviside step function and N is the number of sites. It is easy 
to show that the carrier velocity v (e) is 

« = si = cvm (3) 

The constant C equals (2/\^3j (a/h) (37r 2 )~^ 3 if the system is isotropic in 
three dimensions, a being the lattice constant, given that the charge carriers 
are Fermionic (2 spin values per orbital state). The thermal distribution in (1) 
is classical for the standard reasons that the carriers being injected are so small 
in number, e.g., in time-of-flight observations, that the Fermi distribution is 
well approximated by a Maxwell distribution. Although space-charge limited 
conditions of measurement are common now for some experimental setups such 
as in field-effect-transistor geometry, we will assume here that this assumption 
is valid. If the need arises, it can always be relaxed in subsequent calculations. 

Our present purpose being to get at the simplest consequences of finite B as 
economically as possible, for most of the calculations here, we will first take 
the relaxation time r to have no energy dependence. Such a case appears 
in an early analysis of Erginsoy [15]. Whether the entire set of assumptions 
behind Erginsoy's analysis applies to organic crystals is of no importance in 
the present context. Our interest lies only in exploring, in the first instance, 
the consequences of a constant r in a finite band system. The finite-band 
generalization of (1) is then 

, _ fi /B\ Jo /kBT dxxy/xe- x _ /B\ f B \ $ dx x^xe~ xBlkBT 
* ~ ~ \ SJ tf/ k B T dx y/Ze-* ~ ) \~kjr) lldx^e-^I^T " 



(4) 

Displayed in (4), is the dimensionless ratio // of the mobility to hqC 2 = 
(4a 2 q/3h) (37r 2 ) 2 ^ 3 , and S is the 'scattering energy' H/t. The dimension- 
less mobility // has been expressed in terms of the dimensionless ratios of 
the bandwidth to the scattering energy and to the thermal energy: B/S and 
B/k B T. While both expressions in (4) are exact, they can be used respectively 
in opposite limits most conveniently: 



H 1 = (3/2) (B/S) for B/k B T -> oo 

// = (3/5) (B/S) (B/k B T) for B/k B T -> 0. (5) 

This B to B 2 transition as one goes from large bands to small bands arises 
from the occupation of states in the band. The mobility, which involves the 
thermal average of the square of the velocity, is proportional to Bk B T, the 
product of B and the thermal energy in the large bandwidth case, but to B 2 
in the small bandwidth case, because the temperature is so large with respect 
to B/k B that it makes no contribution: all states in the band are occupied 
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equally. The transition takes on different character if the relaxation time r is 
energy-dependent, the limiting exponents of B being different from 1 and 2 
respectively. 

The two different ways of scaling the energy apparent in (4) are particular 
cases Q = k B T and Q = B of scaling with an arbitrary quantity Q having the 
dimensions of energy: 

' = (!L\ ( Q \ Io B/Q dx ^/xe^f_ 

\Sj \k B T) ffl* dx y/ie-*l**T 1 ] 

This general expression can be used in additional ways such as by taking 
Q = S: 

, _ \ Io /S dx Xy/xe~ xS / kBT 

11 ~ \k^TJ J B/S dx ^e-W ' 

In that case, extreme limits of B/S can be studied — if it is very large, which 
represents the large mean free path case, // is proportional to a ratio of Laplace 
transforms of powers, with S/k B T as the Laplace variable. 

The three representations of the mobility given above are useful in understand- 
ing limiting behavior of carrier transport in extremes of the three respective 
ratios of B, k B T, and S. Generally, one can identify the integrals in (4) with 
incomplete gamma functions 7 defined via 7(6, x) = f$ e~ t t b ~ 1 dt, and write, 
exactly, i.e., for arbitrary relative values of B, k B T, and S, 



/i = 



B\ 7(5/2, B/k B T) 



S) 1 (3/2, B/k B T) 



which can easily be shown to lead to the various respective limits, or rewritten 
in other ways such as in terms only of error functions and exponentials through 
the use of the chain condition 7 (b + 1, x) = 67 (6, x) — x b e~ x and the relation 
7 (l/2,a;) = Vierf( v ^). 



3 Power Laws 



The sharpness of the power laws in the temperature dependence of the ob- 
served fi in organic crystals, whether in the recent experiments on pentacene 
[9] or in experiments on other crystals such as naphthalene reported a couple 
of decades back [12], is impressive. Although Giuggioli et al. [16] have shown 
recently that visually acceptable fits to the data can be produced by a band 
theory addressing the partial range of pentacene data from 20-400 K by com- 
bining acoustic and phonon scattering, it is interesting to ask if a true power 
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law form can be obtained as an analytic limit from the expressions. The follow- 
ing simplified approach produces such as a limit. Assume that the mechanism 
of scattering is such that the relaxation time r depends both on the carrier 
energy e and the temperature T as r = r (e/eo)~ p (T/T )~ r where e , r , 
and T are a characteristic energy, time, and temperature respectively. Such 
expressions arise naturally in many contexts. Thus, p = 1/2 in many semicon- 
ductors [17], and r = 1 for acoustic phonon scattering at high temperatures. 
Equation (1) then yields 

, = (Bto\ ( e \ p 7 (5/2 - p, B/k B T) f T 
* \h)\k B Tj 1 (3/2,B/k B T) \T ) {) 

Whether or not this appears as a power form depends on the particular T 
dependence of the ratio of the 7-functions. For p — 0, one recovers (8) except 
for the explicit T~ r factor. For p—1, the 7-functions cancel and the mobility 
is seen to display a clear temperature law y~( 1+r ). In such the value 

r = 1.7 could conceivably correspond to the observed T~ 2 7 dependence in 
pentacene [9]. It is easy to see that, even if the value of p is not 1, a sharp 
power law may indeed be the direct apparent consequence of the assumed r 
in an appropriate (although limited) range of temperature, provided the ratio 
of the 7-functions in (9) is largely T-independent in that range. 



4 Huang-Rhys Bandwidths and Polaronic Expressions 



If polaronic effects are present in carrier transport, it could be argued [9] that 
some of the physics might be captured by invoking the well-known Huang- 
Rhys dependence 

B = B = B e' G ' 2 c°th(«V2fc B T) (1Q) 

where B is the bare bandwidth, and its exponential reduction occurs as a 
result of strong interaction via the coupling constant G with vibrations of 
frequency Q . For instance, the substitution of (10) in the simple constant-r 
version (8) would yield the normalized mobility 




Equation (11) displays the qualitative trend of the low temperature data, but 
not of the high temperature data, in the acenes. Furthermore, there is always 
a large enough temperature at which the Huang-Rhys factors reduce the po- 
laron bandwidth enough to make the mean free path smaller than a lattice 
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constant. This forces the band calculation to lose its accuracy and, indeed, 
its applicability. Crystal momentum becomes an inadequate quantum number 
in such a case, and a hopping description becomes necessary. In particular, a 
bridging of the present finite-band calculations and of polaronic theories de- 
veloped earlier [4,5] is required. How would one carry out such a combination? 
This important and difficult question, along with application to pentacene ob- 
servations, is being addressed in our ongoing work. Here we sketch a simplified 
answer to the question. 

Polaronic mobility theory for narrow bands has been successfully applied [5] 
more than a decade ago to naphthalene observations [11] on the basis of an 
expression, which, in one of its most simplified forms, may be written as 




where c is a numerical constant and Jo is a modified Bessel function. The 
scattering rate a is essentially the constant 1/r appearing in band expressions 
shown earlier in this Note. The expression comes about from the time inte- 
gration of the velocity autocorrelation function. The velocity autocorrelation 
function is a product of a polaron part which arises from the interaction of the 
charge carrier with vibrations and a part e~ at which arises from interaction 
with static defects [18]. Equation (12) is appropriate to a narrow band. An ex- 
tension of the narrow-band expression to an intermediate-band situation may 
be made by using the calculation given earlier in this Note to obtain the static 
part e~ at . It leads to the replacement of a in (12) by an effective a which is 
dependent on T and B in addition to being proportional to the constant 1/r. 
The result is (with d a numerical constant) 
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5 _B_ 

2' k B T 



3 _B_ 

2' k B T 



sinh (jM=\ 



,2k B Tj . 



The general result (13) combines the polaronic and band-theoretic character 
in a simple way. It reduces to the polaronic form (12) for large T, i.e., for 
k B T » B where B is the reduced bandwidth given by (10). In this limit, the 
reduced band is fully occupied and the ratio of the two 7-functions simplifies 
to (3/5) (B/k B T). In the opposite limit of small temperatures, k B T « B, 
the occupation of the band is controlled by the temperature, the zero-T limit 
of the ratio of the 7-functions is simply 3/2, and (13) simplifies to the band 
result. Indeed, if also k B T << HQ/2, the argument of the J-Bessel function is 
negligible and the Bessel function is simply 1. The mobility then reduces to 
(5) with the zero-T reduced bandwidth B e~ G in place of B in (5). Generally, 
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the mobility, normalized to its T = value, is 



^§- = ^ h iT)lJ 2 f ) (14) 
where \i ah (T) is the right hand side of (11). 

The ad hoc expression (11) obtained by straightforward substitution of the 
bare bandwidth B by the polaronic reduced bandwidth B in (8) differs pre- 
cisely by the introduction in (14) of the Bessel function factor. This agrees 
with the well known fact [4,5] that the excess of that factor over 1 represents 
the hopping contribution of the mobility. We plot (14) in Fig. 1. Parameters 
chosen are arbitrary but illustrate the respective behaviors reported in the lit- 
erature for naphthalene where a flat temperature dependence is observed [11], 
and for pentacene where a turn-over with increasing mobility in the higher 
temperature region dependence has been reported [9]. 



5 Remarks 



Through simplified calculations, we have presented here arguments towards 
an intermediate finite-bandwidth theory which extends previous narrow-band 
treatments. The new results obtained are: finite-band effects stemming from 
full occupation of the band at temperatures large enough so that the thermal 
energy exceeds the bandwidth, B to B 2 transition of the mobility dependence 
on the bandwidth with variation of temperature, an indication of how power 
laws might arise as analytical limits of appropriate scattering mechanisms, the 
specific form of the mobility under the ad hoc Huang-Rhys temperature de- 
pendence of the bandwidths and the relation of such an expression to polaronic 
expressions, and a usable interpolation formula which combines polaronic and 
bare-band characters. 

These calculations should by no means be treated as a complete theory of 
polaronic and band effects. More realistic analyses of the band transport are 
being carried out by Giuggioli et al. [16], taking into account acoustic and op- 
tical phonon interactions which lead to striking new finite-band effects. While 
not totally accurate, the simple truncated free-carrier density of states used 
here has the correct broad (bare) band limit in contrast to Gaussian forms 
[4] and should, therefore, be able to capture the essential physics at that end. 
It is hoped that the present simple calculations will contribute towards the 
construction of a general picture of quasiparticle transport in pure organic 
crystals which treats both the temperature dependence of the mobility, and 
high field effects such as velocity saturation that have been treated recently 
[13]. 



7 



1 



"a 

CD 
N 

"cti 



o 



(b) 



(a) 



Ql , , , 1 

0.5 1 1.5 2 

dimensionless T 

Fig. 1. Dependence of the mobility fi on the temperature T plotted for two (arbi- 
trary) parameter combinations. In units of Ml, the bare bandwidth Bq is 0.5 and the 
coupling constant g is 1.5 in (a) while Bq is 2 and g is 1.8 in (b). The temperature 
T is plotted in units of Ml/2kB and \x is normalized to its value at T = 0. Quali- 
tatively only, (a) and (b) have resemblance to naphthalene and pentacene reported 
data respectively. 
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